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Inspired by the present experimental status of charmed-strange mesons, we perform a systematic study of the 
charmed-strange meson family, in which we calculate the mass spectra of the charmed-strange meson family 
hy taking a screening effect into account in the Godfrey-Isgur model and investigate the corresponding strong 
decays via the quark pair creation model. These phenomenological analyses of charmed-strange mesons not 
only shed light on the features of the observed charmed-strange states, hut also provide important information 
on future experimental search for the missing higher radial and orbital excitations in the charmed-strange meson 
family, which will be valuable task in LHCb, forthcoming Bellell and PANDA. 
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I. INTRODUCTION 

As experiments has largely progressed in the past decade, 
more and more charmed-strange states have been reported [1]. 
Facing the abundant experimental observations, we need to 
provide an answer, as one crucial task, to the question whether 
these states can be identified in the charmed-strange meson 
family, which is not only a valuable research topic relevant 
to the underlying structure of the newly observed charmed- 
strange states, but is also helpful to establish the charmed- 
strange meson family step by step. 

It is a suitable time to give a systematic study of the 
charmed-strange meson family, which includes two main 
topics, i.e., the mass spectrum and strong decay behavior. 
At present, we have abundant experimental information of 
charmed-strange states, which can be combined with the the¬ 
oretical results to carry out the corresponding phenomenolog¬ 
ical study. 

As the first key step of whole study of the charmed-strange 
meson family, the investigation of the mass spectrum of 
charmed-strange mesons should reflect how a charm quark 
interacts with a strange antiquark. Godfrey and Isgur pro¬ 
posed the so-called Godfrey-Isgur (GI) model to describe the 
interaction between q and q quarks inside of mesons [2] some 
thirty years ago. Although the GI model has achieved a great 
success in reproducing/predicting the low lying mesons, there 
exist some difficulties when reproducing the masses of higher 
radial and orbital excitations, which is because the GI model 
is a typical quenched model. A typical example of this de¬ 
fect appears in the low mass puzzle of D*q( 2317) [3-6] and 
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Dji( 2460) [4-7], where the observed masses of D*g(2317) and 
Ds\ (2460) are far lower than the corresponding results calcu¬ 
lated using the GI model. A common feature of higher ex¬ 
citations is that they are near the thresholds of meson pairs, 
which can interact with these higher excitations with the OZI- 
allowed couplings. Hence, it is unsuitable to use the quenched 
GI model to describe the mass spectrum of a higher excitated 
meson and alternatively, we need to adopt an unquenched 
model. Although there have been a couple of works studying 
the heavy-light systems including charmed-strange mesons 
together with their decay modes [8-15], in this work, we 
would like to modify the GI model such that the screening 
effect is introduced to reflect the unquenched pecularity. In 
the following section, we present a detailed introduction of 
the modified GI model. 

In this work, we revisit the mass spectrum of the charmed- 
strange meson to apply the modified GI model , and com¬ 
pare our results with those of the former GI model and ex¬ 
perimental data. We would like to see whether this treatment 
improves the description of the charmed-strange meson spec¬ 
trum to make the modified GI model reliable. We try further 
to obtain information of wave functions of charmed-strange 
mesons, which is important as an input in calculating the de¬ 
cay behavior of the two-body OZI-allowed decay of charmed- 
strange mesons. 

Together with the study of the mass spectrum of charmed- 
strange mesons, it is a valuable information of the property 
of charmed-strange mesons to investigate the decay behav¬ 
ior of charmed-strange mesons. We will adopt the quark 
pair creation (QPC) model [16-22] to calculate the two-body 
OZI-allowed decay of charmed-strange mesons, where the 
corresponding partial and total decay widths are calculated. 
Through this study, we can further test different possible as¬ 
signments to the observed charmed-strange states. In addition, 
we can predict the decay behavior of their partners, which are 
still missing in experiment. This information is important for 
experimentalists to further search for these missing charmed- 
strange mesons, which will be a main task in future experi- 
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ment. 

This work is organized as follows. After Introduction, 
we will briefly review the research status of the observed 
charmed-strange states in Sec. II. In Sec. Ill, the GI model is 
briefly introduced and we give the detailed illustration of how 
to modify the GI model by introducing the screening effect. 
The mass spectrum of charmed-strange mesons together with 
wave fucntions is calculated using the modified GI model. 
Furthermore, the comparison of our results with experimen¬ 
tal data and those from the former GI model is given here. 
With these preparations, in Sec. IV we further study the two- 
body OZI-allowed strong decay behaviors of charmed-strange 
mesons through the QPC model, where the QPC model is 
briefly introduced. Next, we perform the phenomenological 
analysis by combining our results with the experimental data. 
The paper ends with a summary in Sec. V. 

II. CONCISE REVIEW OE THE OBSERVED 
CHARMED-STRANGE STATES 

In this section, we briefly review the experimental and theo¬ 
retical status on charmed-strange mesons. First, in Table I we 
collect the experimental information of the observed charmed- 
strange states, which include resonance parameters and the 
corresponding experiments. Since D* and D** have beeen es¬ 
tablished to be IS cs mesons, our review mainly focuses on 
possible candidates of higher radial and orbital excitations in 
the charmed-strange meson family. 

A. 0,1 (2536) and D;2(2573) 


dominantly occurs via 5-wave with F(Dji(2536)^ ^ 
£>*^7i:«)s_wave/r(D,i(2536)+ ^ = (0.72 + 0.05 + 

0.01). The mass value and narrow width of Dji(2536) are con¬ 
sistent with the theoretical expectation of the charmed strange 
meson as 7^ = N [36]. 

In 1994, the CLEG Collaboration observed a charmed- 
strange meson 0*2(2573) in the invariant mass spec¬ 

trum [37], where the measured resonance parameters are M = 
(2573+|-^+0.8+0.5) MeV and F = (16+5+3) MeV. In addition, 
an upper limit of the following ratio was given 

S(0*A2573)+ ^ D*^K+) 

-- - < 0.33, (2) 

S(0*2(2573)+ ^ DPK+) 

which is from the search for the decay mode 0*2(2573)+ —> 
O*o^+ pyj present, 0*2(2573) is a good candidate of a 
charmed-strange meson with \^P 2 since its mass is consistent 
with the theoretical prediction [36]. The 0*2(2573) had sub¬ 
sequently been confirmed by some other collaborations in the 
0°^r+ invariant mass spectrum [29, 38-41]. 


B. 0]„(2317) and 0,1(2460) 

In 2003, the BaBar Collaboration observed a new charmed- 
strange state near 2.32 MeV in the O+itt** invariant mass dis¬ 
tribution of the B decay, which is named as 0*^(2317) [3]. 
Later, the CLEG Collaboration confirmed this state in the 
0 + 7 r° channel and reported another narrow charmed-strange 
state 0,1 (2460) [4]. In addition, the ratio 


Before 2013, there are only two good candidates for the 
IP states in the charmed-strange meson family, which are 
0,1 (2536) and 0^2(2573). 

As a charmed-strange meson with 7^ = 1+, 0,i(2536) was 
first observed in 1987 by analyzing the 0*7 invariant mass 
spectrum of the vN scattering process [23], where the mea¬ 
sured mass is (2535 ± 28) MeV. Later, the ARGUS Collabora¬ 
tion observed this state in the D*^K^ final state, where its mass 
and width are M = (2536 + 0.6 ± 2.0) MeV and F < 4.6 MeV 
[24]1, respectively. Since there does not exist the 0,i(2536) 
signal in the O+^f*^ invariant mass spectrum, 0,i(2536) has an 
unnatural spin-parity [24]. In 1993, the CLEG Collaboration 
measured the ratio of F(0,i(2536) —> 0 * 7 ) to F(0,i(2536) ^ 
D*K), which is [25] 


F(0,i(2536) ^ 0 : 7 ) 
F(0,i(2536) ^ D*K) 


( 1 ) 


The 0,1 (2536) has been confirmed by other groups in different 
channels [26-34]. The BaBar Collaboration reported this state 
in the O+tt+tt^ invariant mass spectrum [ 6 ]. The Belle Col¬ 
laboration observed 0,i(2536) —> D^n K^ decay mode [35], 
where the ratio S(0+j(2536) —> 0+7r"/r+)/S(0+j(2536) ^ 
0*+ir+) = (3.27 + 0.18 + 0.37)% was obtained. In 
addition, the measurement of the angular distribution of 
the 0,1 (2536)+ —» D*^K^ indicates that 0,i(2536)+ —> 


S(0*o(2317)+ ^ 0:+7) 
S(0;o(2317)+ ^ 0 + 7 r«) 


< 0.059 


(3) 


was obtained in Ref. [4] while the Belle and BaBar collab¬ 
orations gave the upper bound of this ratio as 0.18 and 0.16, 
respectively [5, 6 ]. The O,i(2460) was also confirmed by the 
Belle and BaBar collaborations [5-7]. 

If assigning 0*^(2317) and O,i(2460) as charmed-strange 
mesons with quantum numbers 7^ = 0 + and 7^ = 1 +, re¬ 
spectively, there exists the low mass puzzle for 0*^(2317) 
and 0,1 (2460), i.e., the theoretical masses of charmed-strange 
mesons with 0+ and 1+ are 2.48 GeV and 2.53 GeV, respec¬ 
tively, predicted by the old but successful model proposed 
in Refs. [2, 36], which values are far larger than the cor¬ 
responding experimental data. These peculiarities have also 
stimulated theorists’ extensive interest in exploring their inner 
structures and the exotic state explanations to 0 *q( 2317) and 
0,1 (2460) were especially proposed in Refs. [43-45]. 

Since 0 *q( 2317) and O,i(2460) are near and below the 
thresholds of DK and D*K, respectively, the coupled-channel 
effect, which is an important non-perturbative QCD effect, 
should be considered in understanding the low mass puzzle. 
This means that 0*^(2317) and O,i(2460) are still categorized 
as the conventional charmed-strange meson family [46-48]. 
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TABLE I: Experimental information of the observed charmed-strange states. 


State 

Mass (MeV) [1] 

Width (MeV) [1] 

1” observation 

Observed decay modes 

D, 

1968.49 ± 0.33 




d: 

2112.3 + 0.5 

< 1.9 



D‘J23n) 

2317.8 + 0.6 

< 3.8 

BaBar [3] 

D)7r° [3] 

£>.u(2460) 

2459.6 ± 0.6 

< 3.5 

CLEO [4] 

D;+7r° [4] 

T>.u(2536) 

2535.12 + 0.13 

0.92 + 0.05 

ITEP, SERF [23] 

D^y [23] 

£ 1 : 2 ( 2573 ) 

2571.9 + 0.8 

16!J + 3 [37] 

CLEO [37] 

D°K* [37] 

L»:/2632)“ 

2632.5 + 1.7 ±5.0 [40] 

< 17 [40] 

SELEX [40] 

D°K+ [40] 

£':,(2700) 

2688 ±4 + 3 [39] 

112 + 7 + 36 [39] 

BaBar [39] 

DK [39] 


2708 + 9];]* [49] 

108 + 23!“ [49] 

Belle [49] 

D°K* [49] 


2710±2];]2 [50] 

149 + 7!^^ [50] 

BaBar [50] 

£)<*>X[50] 


2709.2+ 1.9 + 4.5 [51] 

115.8 + 7.3+12.1 [51] 

LHCb [51] 

DK [51] 

DIj(2860) 

2856.6+ 1.5 + 5.0 [39] 

47 + 7+10 [39] 

BaBar [39] 

DK [39] 


2862 + 2)1 [50] 

48 + 3 + 6 [50] 

BaBar [50] 

[50] 


2866.1 + 1.0+ 6.3 [51] 

69.9 + 3.2 + 6.6 [51] 

LHCb [51] 

DK [51] 

£>:3(2860) 

2860.5 + 2.6 + 2.5 + 6.0 [62, 63] 

53 + 7 + 4 + 6 [62, 63] 

LHCb [62, 63] 

rfK- [62, 63] 

£):i(2860) 

2859 + 12 + 6 + 23 [62, 63] 

159 + 23 + 27 + 72 [62, 63] 

LHCb [62, 63] 

D'^K- [62, 63] 

O.y(3040) 

3044 + 8];f [50] 

239 + 35!“ [50] 

BaBar [50] 

D-K [50] 


“Since D*j{2632) was only reported by SELEX [40] and not confirmed by 
other experiments [42], we do not include this state in our review of this work. 


C. £>:i(2700) 


In 2006, the BaBar Collaboration observed a broad struc¬ 
ture in the DK invariant mass spectrum, which was named as 
Z)*j(2700) [39] and the mass and width are M - (2688+4 + 3) 
MeV and F = (112 + 7 + 36) MeV, respectively. Later, 
Z)*j(2700) was conhrmed by the Belle Collaboration in the 
DK invariant mass spectrum of B ^ D'^{D^K^} process [49]. 
The angular momentum and parity of D*j(2700) are deter¬ 
mined to be 7 = 1 and P = - by the helicity angle dis¬ 
tribution and by its decay to two pseudoscalar mesons, re¬ 
spectively. The Babar Collaboration reported the D*K de¬ 
cay mode of Dji(2700) in Ref. [50], and obtained a ratio of 
S(D;j(2710) ^ D*K) toS(D;;(2710) ^ DK) as [50] 


S(D;j(2710) ^ D*K) 
S(D;j(2710) ^ DK) 


0.91+0.13 + 0.12. 


(4) 


In 2012, the LHCb Collaboration also observed D*j(2700) in 
the DK mass spectrum [51]. 

As a vector charmed-strange state, the measured mass of 
Z)*j(2700) is close to the prediction of the 2^5 1 charmed- 
strange meson [2]. The strong decay behavior of a cs state 
of 2^5 1 is investigated using the QPC model in Ref. [52]. 
The ratio S(D*j(2700) ^ D*K)IS(D*J2700) DK) evalu¬ 
ated by the effective Lagrangian approach, however, favors the 
2*5o assignment for D*j(2700) [53]. Besides the decay be¬ 
havior, the production of D*j(2700) from the B meson decay 
was calculated by a naive factorization method, which shows 
that D*j(2700) could be explained as the first radial excitation 
ofD*,(2112) [54]. 


In Refs. [55, 56], D*j(2700) was assigned as a mix¬ 
ing of 2^5 1 and l^Di cs state, and this assignment can be 
supported by the study of the decay behavior obtained by 
the QPC model. The obtained ratio of S(D*j(2710) —> 
D*/r)/S(D*j(2710) —> DK) is also consistent with the exper¬ 
imental measurement [56] from BaBar [50]. In addition, eval¬ 
uation by the constituent quark model [57] and the Regge Phe¬ 
nomenology [60] also support the assignment of D*j(2710) as 
a mixing of the 2^S \ and l^Di cs states. Other than a stan¬ 
dard interpretation of D*j(2700) as a cs state, a molecular state 
explanation was proposed in Ref. [61], which is based on a 
potential model. 


D. £);/2860), (2860) and (2860) 


The D*^(2860) was first discovered by the BaBar Col¬ 
laboration in the DK invariant mass spectrum of the inclu¬ 
sive process —» DKX [39]. Its resonance parame¬ 

ters are reported d& M — (2856.6 + 1.5 + 5.0) MeV and 
r = (48 + 7+10) MeV. This state was conhrmed in the 
D*K mode by the BaBar Collaboration [50] again and the ra¬ 
tio S(DJ^(2860)+ ^ D*/:)/S(D;^(2860)+ ^ DK) was mea¬ 
sured as 


S(D;^(2860)+ ^ D*K) 
S(D;^(2860)+ ^ DK) 


1.10 + 0.15 + 0.19. 


(5) 


Very recently, the LHCb Collaboration reported their new 
measurement of the structure around 2.86 GeV in the D^K^ 
invariant mass distribution of —> D^K tt^ [62, 63]. The 
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amplitude analysis of this decay indicates that the structure 
at - 2.86 GeV contains both - 1“ and 3“ com¬ 

ponents corresponding to D*j(2860) and D*^(2S6Q), respec¬ 
tively, where the resonance parameters of D*j(2860) and 
D* 3 ( 2860 ) are 


4Ly)*j(2860) 

= (2859+ 12 + 6 + 23 )MeV, 

( 6 ) 

r£);,( 2860 ) 

= (159 + 23 + 27 + 72)MeV, 

(7) 

MdI^(2S60) 

= (2860.5 + 2.6 + 2.5 + 6.0 ) MeV, 

( 8 ) 

rD; 3 ( 2860 ) 

= (53+7 + 4 + 6 )MeV. 

(9) 


Before the measurement by the LHCb Collaboration, the 
properties of D*^(2860) have been widely discussed. The pos¬ 
sibility of D*y(2860) as the first radial excitation of D*q( 2317) 
has been ruled out due to the observation of D*y(2860) ^ 
D*K [50]. According to the calculation by the QPC model 
[52], the Regge phenomenology [60], chiral quark model 
[59], and the flux tube model [64], D*y(2860) can be as¬ 
signed to a charmed strange meson. However, dif¬ 

ferent approaches give different ratios of S(D*y(2860) ^ 
D*K)I!B(D*j(2S60) —> DK). For example, this ratio was esti¬ 
mated to be 0.36 by the effective Lagrangian method [53] and 
to be 0.8 by the QPC model [56]. At present calculation, the 
= 3“ assignment to D*y(2860) is still possible. 

In Ref. [56], the 2S-1D mixing was proposed to ex¬ 
plain D*^(2860), where D*j(2860) and D*j(2710) are treated 
as a mixture of 2^5 1 and l^Di states in the charmed- 
strange meson family. With a proper mixing angle, the ra¬ 
tios of S(D*j(2860) D*K)IS(DIj{2860) DK) and 

!B{Dlj(2700) D*K)IS{Dlj(2700) DK) can be well ex¬ 

plained simultaneously. In Ref. [57] the authors proposed a 
two-state scenario for D*^(2860), one resonance is likely to be 
the l^D^ and the other resonance seems to be the higher mix¬ 
ing state of 1 ^D 2 - I'Di- Although the - 0^ assignment 
to D*j{286Q) has been ruled out by the experimental measure¬ 
ment of D*j(2860) —» D*K, the authors of Ref. [65] indicated 
that there exist two resonances around 2.86 GeV with - 0^ 
and - 2^. In the DK invariant mass spectrum, the struc¬ 
ture near 2.86 GeV should contain these two resonances while 
in the D*K invariant mass spectrum, only one resonance of 
= 2 ^ should be included. 

Following the measurement by the LHCb Collaboration 
[62, 63], the decay behaviors of D*j(2860) and D*3(2860) 
were evaluated by the QPC model [ 66 , 68 ], which show that 
these states can be good candidates of the ID states in the 
charmed-strange meson family. By using the QCD Sum Rule, 
the masses of ID charmed-strange mesons were calculated 
in Ref. [69], which also supports the explanation of the ID 
charmed-strange mesons. In addition, another study of the 
decay behaviors of theses states were performed in Ref. [67], 
where the effective Lagrangian approach was adopted. 

Before closing the review of D*y(2860), Z)*j(2860) and 
D* 3 ( 2860 ), we need to comment on the measurement of the 
ratio in Eq. (5). Since the new measurement by LHCb [62,63] 
is given, thid ratio must be changed according to which state 
is assigned to D*j in both denomenator and numerator. Thus, 
we do not suggest to adopt the old data of this ratio in Eq. (5) 
when checking the theoretical result. We also expect new 


measurement of this ratio when considering the LHCb results 
[62, 63], which will be helpful to pin down the different ex¬ 
planations. 


E. Z),y(3040) 

Besides confirming D*j(2700) and D*^(2860) in the D*K 
invariant mass distribution, the BaBar Collaboration also ob¬ 
served a new broad structure with mass M - (3044 + 8 ]; 5 °) 
MeV and width F = (239 + 35 ^^ 42 ) MeV. The negative result 
of a decay Djy(3040) —> DK suggests unnatural parity for 
£>.y(3040). 

The observed mass of Dsy(3040) and its unnatural par¬ 
ity are consistent with the quark model prediction for 2 ^Pi 
charmed strange meson [ 2 ], which is the first radial excitation 
of Dii(2460). The calculations in the QPC model also support 
that Djy(3040) can be categorized as 1^ state in a (0^, 1^) spin 
doublet [70] of the heavy quark symmetry. In addition, the 
calculations in the flux tube model [64], the constitute quark 
model [57, 58] and the effective approach [71] also indicate 
that a possibility of a Dsy(3040) as a 1'^ charmed-strange me¬ 
son can not be ruled out. References [72, 73] calculated the 
decay width of Djy(3040) as n(J^) - 3(1^) or 4(1^) state, 
which is rather large but still compatible with the experimen¬ 
tal data. Besides 7^ = 1^ assignment to Djy(3040), the possi¬ 
bility that D5y(3040) as a mixture of l^Dy and 1 'Dy charmed- 
strange mesons was discussed with an effective Lagrangian 
approach [71]. 

As can be seen in the above review of status of the observed 
charmed-strange states, we notice that more and more can¬ 
didates of higher radial and orbital charmed-strange mesons 
were reported in the past decade. It is a good opportunity 
to carry out the systematic study of charmed-strange mesons 
now, which will deepen our understanding of the charmed- 
strange meson family and will provide more abundant infor¬ 
mation for experimentalists to further search for higher radial 
and orbital charmed-strange mesons in future experiment. 

In the following sections, we start to investigate charmed- 
strange mesons on their mass spectrum and two-body OZI- 
allowed decay behaviors. 


III. MASS SPECTRA 

In this work, we employ the modified relativistic quark 
model to calculate the mass spectra and wave functions of 
charmed-strange mesons because owing to the peculiarity of 
charmed-strange mesons, relativistic effects and unquenched 
effects cannot be ignored. In 1985, Godfrey and Isgur pro¬ 
posed the GI model to describe the meson spectra with a great 
success, especially for the low lying mesons [2] and later God¬ 
frey and Kokoski studied E-wave heavy-light systems [36]. 
Since a coupled-channel effect becomes more important for 
higher radial and orbital excitations, we need to include this 
effect in calculating the mass spectrum, which motivates us to 
modify the GI model by introducing the screened potential. 
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which is partly equivalent description to the coupled-channel 
effect [74], 

In the following, we first give a brief review of the GI model 
and then present how to introduce the screened potential. 


A. Brief review of the Godfrey-Isgur model 

The interaction between quark and antiquark of the GI 
model [2] is described by the Hamiltonian 

^ = (/ + + (/ + '" 2 )*^^ + (p, r), ( 10 ) 

where yeff(p, r) = + ffso effective potential 

of the qq interaction. Veff(p,r) contains two main ingredients. 
The first one is a short-distance ® 7 ^ interaction of one- 
gluon-exchange and the second is a long-distance 1 ® 1 linear 
confining interaction which is at first employed byt the Cor¬ 
nell group and is suggested by the lattice QCD. This effective 
potential can be obtained by on-shell qq scattering amplitudes 
in the center-of-mass (CM) frame [2]. 

In the nonrelativistic limit, Veff(p, r) is tansformed into the 
standard nonrelativistic potential Veff(r). 

V^s(r) = H- + H^° ( 11 ) 


with 


//“"f ^c + br+ ■ F 2 , (12) 

r 

where is the spin-independent potential and contains 
a constant term, a linear confining potential and a one-gluon 
exchange potential. The subscripts 1 and 2 denote quark and 
antiquark, respectively. The second term of r.h.s. ofEq. (11) 
is the color-hyperfine interaction, i.e.. 


//hyp - 


a sir) 


mym2 


Stt 




1 /3S'i ■ rS '2 ■ r 


-Si ■ S 2 ) 


Ft ■ F 2 . 


(13) 


The third term of r.h.s. of Eq. (11) is the spin-orbit interaction, 
//SO _ ^SO(cm) _j_ //SO(tp) (14) 

Here, //SO(cm) jj^g color-magnetic term and is the 

Thomas-precession term, i.e.. 


^SO(cm) _ _ °^sir) ^ 


^SO(tp) _ 



m\ m 2 /\mi m 2 ) 


2r dr mj) 


■ L(Fi ■ F 2 ), (15) 


L. (16) 


quark and antiquark. F is related to the Gell-Mann matrix by 
Fi = Ai/2 and F 2 = Eor a meson, {F\ ■ F 2 ) - -A 13. 

The GI model constructed by Godfrey and Isgur is a rel¬ 
ativized quark model, where relativistic effects are embeded 
into the model mainly in two ways. 

Eirstly, a smearing function pnir — r') is introduced to 
incorporate effects of an internal motion inside a hadron and 
nonlocality of interactions between quark and antiquark. A 
general form is given by 


fir) 




'piiir - r')f(r') 


(17) 


with 


Pi2(r - r') = ^exp 


{-(Aiiir - r’f ), 


(18) 


^12 “ 


1 1 / 4mim2 

2 2\(miH-OT2)^ 


-I- s 


2 j 2m\m2 
m\ -I- m 2 


where cro — 1 .80 GeV and s — 1.55 are the universal param¬ 
eters in the GI model, mj and m 2 are masses of quark and 
anti quark, respectively. If nii = m 2 and both of them become 
large, cri 2 also becomes large andpi 2 (r—r') ^ d^ir—r') and 
hence this smearing is effective for a heavy quarkonium be¬ 
cause the second term becomes (smif. If m 2 » mi, cri 2 gets 
its minimum and suitable for describing a heavy-light system 
like the charmed-strange mesons. 

Secondly, a general formula of the potential should depend 
on the CM momentum of the interacting quarks. This effect is 
taken into account by introducing momentum-dependent fac¬ 
tors in the interactions and the factors will go to unity in the 
nonrelativistic limit. In a semiquantitatively relativistic treat¬ 
ment, the smeared Coulomb term G(r) and the smeared hy- 
perfine interactions V, should be modified according to. 


Gir) 


Viir) 

mim2 


1 -H 


E 1 E 2 

, 1 / 2 +fi 


1/2 


<5(0 1 


I mim 2 \ Vijr) / otiOT 2 \ 

\£'i£' 2/ miOT2\£'i£'2/ 


E 1 E 2 ) 

1/2+ei 


1/2 


(19) 


where Ei and £2 are the energes of the quark and antiquark in 
the meson, e, ’s are expected to be small numbers for a different 
type of hyperfine interactions, i.e., the contract, tensor, vector 
spin-orbit and scalar spin-orbit potentials, into which Yfr) are 
classified. The particular values of e, and other parameters are 
given in Table II of Ref. [2]. 

Diagonalizing the Hamiltonian given by Eq. (10) by a sim¬ 
ple harmonic oscillator (SHO) basis and using a variational 
method, we get the mass spectrum and wave function of a me¬ 
son. More details of the GI model can be found in Appendices 
of Ref. [2]. 


B. Modified GI model with the screening effect 


In the above expressions, S 1 /S 2 denotes the spin of the Although the GI model has achieved a great success in de¬ 
quark/antiquark and L is the orbital momentum between scribing the meson spectrum [2], there still exists discrepancy 
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between the predictions given by the GI model and recent 
new experimental observations. A typical example is that the 
masses of the observed Dso{22>ll) [3-6], Dji(2460) [4-7] and 
X(3872) [76] deviate from the corresponding values expected 
by the GI model [2]. Later, theorists realized that these dis¬ 
crepancies are partly caused by coupled channel effects [46- 
48, 77]. Furthermore, another remedy can be made to ad¬ 
just masses by screening the color charges at distances greater 
than about 1 fm [78] which spontaneously creates light quark- 
antiquark pairs. The screening effect has been conhrmed by 
unquenched Lattice QCD and some holographic models [79- 
82]. 

There were several progresses on the study of the meson 
mass spectrum by considering the screening effect [83, 84]. In 
Ref. [83], the authors adopted the screened potential to com¬ 
pute the charmonium spectrum [74]. Mezzoir et al. in [84] 
carried out the investigation of highly excited light unflavored 
mesons by flattening the linear potential br above a certain 
saturation distance r^. However, the study of a screening ef¬ 
fect for a heavy-light meson system is still absent at present, 
which is the reason why we introduce the screening effect into 
the GI model in this work. 

In order to take into account the screening effect, we need 
to make a replacement in Eq. (12) [85, 86] as 


br ^ = 




( 20 ) 


where V^^’^(r) behaves like a linear potential br at short dis¬ 
tances and approaches to ^ at long distances. Then, we further 
modify y^‘^'^(r) as the way given in Eq. (17), 

d^r'pnir — r')— -. (21) 

By inserting Eq. (18) into the above expression, the concrete 
expression for y®‘^' (r) is given by 


y^“(r) = - 

pr 


r -I- e 


Jn Jo 2 


2cr2 \ 




p - 2rcr^ j 

^ 1 

2cr2 1 



—— I e ^ dx - - 


=r 


( 22 ) 


with 
1 d 


r or 


fin 


r + 2r^a-^p-p 
2cr^r 


p+lra^ 


e -"dx - ^ 1 + 


ur-2rcrp+p l 


2cr^r 




\ 


Lf- 

Jn Jo 




e ^ dx 
(23) 



r(fm) 


FIG. 1: (color online). The r dependence of V(r) = br and 

V(r) = f cPr'p[ 2 (r — r')br'. Here, we take two values p = 0.02 
GeV and p = 0.05 GeV to show y®“(r). Other parameters involved 
in y*“(r), V(r) = br and V{r) are m., = 1.628 GeV, m, = 0.419 GeV, 
(To = 1.80 GeV and i = 1.55 in the expression of (T| 2 . See TABLE 
II in [2]. 


In Eig. 1, we compare the r dependence of our yscr(^)^ usual 
linear potential y(r) = br and y(r) = J d^r'p\ 2 (r — r')br'^. 
While behaviors of W“(r) and V(r) are similar to each other 
at small distances, there exists difference between y(r) = br 
and smeared Vir) reflecting the relativistic effect. There is 
obvious difference between y^‘^'^(r) and y(r) at large distances 
due to the screening effect. Comparing a heavy-light system 
with a heavy quarkonium, a velocity of a light quark in the 
heavy light meson is much larger than that of a heavy quark 
in the quarkonium and in addition, a radius of the heavy light 
meson is also larger than the heavy quarkonium. These facts 
indicate that we need to consider the relativistic and screening 
effects when studying the mass spectra of higher radial and 
orbital excitations of charmed-strange mesons. Eurthermore, 
adjusting p is to control the power of screening effect. 


C. Numerical results 

In the following, we present numerical results of the mass 
spectrum of the charmed-strange meson family. In Table II, 
the masses of the charmed-strange mesons are listed using the 
GI model [2] and our modified GI model, where we take sev¬ 
eral p values in W‘^'(r) to show p dependency of our model. 
Apart from p, other parameters appearing in our calculation 
are taken from Ref. [2], with which we can, of course, repro¬ 
duce the masses obtained in Refs. [2, 75]. In Table II, we also 


* The concrete expression for y(r) is given in Ref. [2] as 


Jncrr 


+ 1 + 



e dx 


where V(r) is given in the footnote in this page and cr - (Ti 2 
given by Eq. (18). 


V(r) = br 























7 


TABLE II: Mass spectrum of the charmed-strange meson family obtained hy our modified GI model and the comparison with those calculated 
by the GI model. Here, we take yu = 0.01,0.02,0.03, and 0.04 GeV to show the results with the modified GI model. The values in brackets are 
the obtained R= Ijp and n denotes the radial quantum number. We emphasize that we do not consider the mixing among states with the same 
quantum number when presenting the results in this Table, ji is in units of GeV, while R is units of GeV^'. 


GI model [2, 75] 

n = 1 « = 2 

n = 0.01 

n = 

^ = 0.02 

1 

H = 0.03 

Modified GI model 

^ = 0.04 p = 0.01 

n = 

11 = 0.02 

2 

fi - 0.03 

p = 0.04 

n'So 

1979(1.41) 

2673(2.00) 

1971 

1967(1.41) 

1963 

1960 

2661 

2646(2.08) 

2632 

2618 

n^S 1 

2129(1.69) 

2732(2.08) 

2120 

2115(1.71) 

2111 

2106 

2719 

2704(2.17) 

2688 

2673 

n'Pi 

2550(1.92) 

3024(2.17) 

2541 

2531(1.96) 

2522 

2512 

3001 

2979(2.27) 

2957 

2935 

n^Po 

2484(1.75) 

3005(2.13) 

2473 

2463(1.79) 

2454 

2444 

2982 

2960(2.22) 

2937 

2914 

n^P, 

2552(1.89) 

3033(2.22) 

2542 

2532(1.96) 

2522 

2512 

3010 

2988(2.27) 

2965 

2942 

n^P2 

2592(2.08) 

3048(2.27) 

2581 

2571(2.17) 

2561 

2551 

3026 

3004(2.38) 

2981 

2959 


2910(2.22) 

3307(2.38) 

2893 

2877(2.27) 

2861 

2844 

3277 

3247(2.50) 

3216 

3186 

n^Di 

2899(2.08) 

3306(2.27) 

2882 

2865(2.13) 

2848 

2831 

3275 

3244(2.44) 

3213 

3182 

n^D2 

2916(2.17) 

3313(2.38) 

2899 

2882(2.27) 

2865 

2848 

3283 

3252(2.50) 

3221 

3190 


2917(2.33) 

3311(2.44) 

2900 

2883(2.38) 

2867 

2850 

3281 

3251(2.56) 

3221 

3190 


3199(2.38) 

- 

3175 

3151(2.50) 

3127 

3102 

- 

- 

- 

- 

n^F2 

3208(2.27) 

- 

3183 

3159(2.38) 

3134 

3109 

- 

- 

- 

- 


3205(2.33) 

- 

3181 

3157(2.44) 

3132 

3107 

- 

- 

- 

- 

n^F^ 

3190(2.44) 

- 

3167 

3143(2.56) 

3120 

3096 

- 

- 

- 

- 


TABLE III: Comparison of the experimental data and theoretical results. Here, we also list thex~ values for different models. The notation 
is introduced to express mixing states either of 'L/. or ^L/,. 



n 

Experimental values [1] 

GI model [2, 75] Modified GI model" 

D. 

1 'So 

1968.49 ± 0.33 

1979 

1967 

d: 

1^5, 

2112.3 ±0.5 

2129 

2115 

H,i(2536) 

ITi 

2535.12 ±0.13 

2556 

2534 

0 : 2 ( 2573 ) 

i^Pi 

2571.9 ±0.8 

2592 

2571 

0:^(2700) 

2 351 

2709 ± 4 

2732 

2704 

£>:i(2860) 

I'Oi 

2859± 12±6±23 [62, 63] 

2899 

2865 

£>:3(2860) 

I'Os 

2860.5 ± 2.6 ± 2.5 ± 6.0 [62, 63] 

2917 

2883 

Z),,(3040) 

2P, 

3044 ± 8!f [50] 

3038 

2992 


- 

- 

7165 

36 


“The results listed in the last column are calculated via the modified GI 
model with fj = 0.02 GeV. 

give the R - lf/3 values corresponding to charmed-strange 
mesons, where the R value can be obtained by the relation 

f - f 0(p)2/^3p (24) 

The r.h.s of Eq. (24) is the root mean square momentum, 
which can be directly calculated through the GI model or the 


modified GI model. The I.h.s of Eq. (24) is a definition of the 
root mean square momentum when adopting the SHO wave 
function. Here, in the momentum space the SHO wave func¬ 
tion is expressed as 


« (p) = P!^°(p)rLM,(n,), (25) 
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FIG. 2: (color online). Mass spectrum of charmed-strange mesons (in units of MeV). Here, the blue lines stand for the results obtained by the 
GI model [2, 75], while the red lines are those by the modified GI model with p = 0.02 GeV. The black rectangles denote the experimental 
data taken from PDG [1]. The corresponding quantum numbers are listed in the abscissa. In addition, when there exists a mixture of 

and n^Li states, we use another notation L^. 


with the radial wave function 


defined as 


RT(p) 


i-m-if 


2n\ 


^03/2 y^T(n + L + 3l2)\l3 


P\ 

‘ e 2 /j2 


(26) 


where 10^) is an associated Laguerre polynomial 

with the oscillator parameter p. In the configuration space 
the SHO wave function is given by 

= Rir(r)YLM0n,), (27) 


where the radial wave function is defined as 




= «3/2 


2 «! 


I r(n -I- L -H 3/2) 


(Pr0e L0^l^(p^0).{22,) 


The p is the parameter appearing in the SHO radial wave func¬ 
tion given in Eq. (26), which is determined by the above pro¬ 
cedure. 

In Table III, we further compare the calculated results with 
the experimental data. Illustrating further a suitable p value 
introduced in the modihed GI model, we also list the val¬ 
ues for the GI model and the modihed GI model, where is 




2 



-^Th(0 - -^Exp(0 

Error(0 


(29) 


where ^Th(i) and 2^Exp(i) are theoretical and experimental 
values. Error(0 is the experimental error listed in Table III. 
In Table III, the values of x^ are calculated without the con¬ 
tributions from Djo( 2317) and Dji(2460). Comparing the;^^^ 
values for different cases, we hnd that the modihed GI model 
can well describe the experiments when p = 0.02 GeV and 
that the description of the observed charmed-strange spectrum 
can be obvisously improved by the modihed GI model with 
the screening effect. One can clearly see the trend of decreas¬ 
ing masses with increasing p from Table III. 

In Eig. 2, we further list the results obtained by the modihed 
GI model (p - 0.02 GeV) and make a comparison with the 
GI model and experimental data. As for Ds, Z)*, Dji(2536), 
D*2(2573), Z)*j( 2700) and the newly observed D*j(2860), the 
mass differences between the experimental data and our theo¬ 
retical results are less than 10 MeV. While the theoretical mass 
of 0*^(2860) is about 20 MeV lower than the experimental 
value, we notice that the experimental masses of D*q( 2317) 
and Dji( 2460) can not be reproduced by the modihed GI 
model. This discrepancy is caused by the near-threshold ef- 
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feet which is ignored in the screening effect [74]. In general, 
the comparison of 10 experimental data with our calculated re¬ 
sults indicates that the modified GI model is more suitable to 
describe the experimental data, especially to higher charmed- 
strange mesons. 

We should emphasize that the n'Pi -n^Pi and n'Z )2 -n^Z )2 
mixtures with n - 1,2 are included in the corresponding cal¬ 
culations listed in Table III and Fig. 2, which is different form 
the situation in Table II since we need to compare theoretical 
results with the experimental data. 


A as an example, we illustrate the definition of a mock state, 
i.e., 

\A{n^^^^LjM,){VA)) 

- 2 {LMl-SMs\JMj)x^sm. 

Ms,M l 

X<P^CJ^ Jdpidp 2 S^(PA - Pi - P 2 ) 
X'i'«LMi.(Pl-P2)l?l(Pl)?2(P2)), (32) 


IV. TWO-BODY OZI-ALLOWED STRONG DECAYS 

In addition to the mass spectrum, the decay properties are 
also crucial features of mesons. The quarks pair creation 
(QPC) model is successful to calculate Okubo-Zweig-Iizuka 
(OZI) allowed strong decays of mesons. Here, we firstly give 
a brief introduction to the QPC model. 


A. Brief introduction to the QPC model 

The QPC model was proposed by Micu [16] and developed 
by the Orsay group [17-22]. It assumes that meson decay oc¬ 
curs through a flavor-singlet and color-singlet quark-antiquark 
pair created from the vacuum. The transition operator T can 
be expressed as 


where Xs m ’ flavor and color wave 

functions, respectively. P 2 ) the spacial wave 

function of a meson A in the momentum space, which can be 
obtained by the modified GI model. The calculated amplitude 
]\^M]aMjbM]c converted into the partial wave amplitude 

M-tt- via the Jacob-Wick formula [89], i.e.. 


M'^A BC) 


V2L H- 1 1 

—— 2 ] mJMjA\JAMjA) 

^ MjbMjC 

x{JbMjb; JcMjc\JMja)M^^^^^‘^^^. 

(33) 


Finally, the decay width can be expressed as 


F = 



J,L 


(34) 


T - -3y 2^(lm; 1 - m|00) j t/p3t/p45^(P3 + P 4 ) 

X (30) 

where p 3 and p 4 denote the momenta of quark and anti¬ 
quark created from the vacuum, respectively. 7 is the pa¬ 
rameter which describes the strength of the creation of quark- 
antiquark pair. By comparing the experimental widths with 
theoretical ones of 16 famous decay channels, y = 8.7 is ob¬ 
tained for the uuidd pair creation [87]. For the ss pair cre¬ 
ation, we take y - 8.7/ V3 [17]. - (uu + dd + ss)/ ^/3 is 

the flavor function while = dij/ V3 is the color func¬ 

tion, where i and j are the color indices, iftmip) - Ipl^Tfm(p) 
denotes the solid harmonic polynomial. Thus, J/im(®^ 2 ^) in¬ 
dicates that the angular momentum of the quark-antiquark pair 
is L = 1. The means that the total spin angular momen¬ 
tum of the quark-antiquark parr is 5 = 1. Finally, the quantum 
number of the quark-antiquark pair is contained as 7^'' = 0 ^^ 
through coupling of the angular momentum with the spin an¬ 
gular momentum. 

The transition matrix of meson A decaying into mesons B 
and C in the A rest frame is defined as 

{BC\T\A) = 5^(pb -h pc)yV(4rMM,«Mjc^ (31) 

Pb and pc are the momenta of mesons B and C, respectively. 
|A), \B} and |C) denote the mock states [ 88 ]. Taking a meson 


We would like to emphasize the improvement of the present 
work compared with the former works in Refs. [52, 66 , 70]. 
In Refs. [52, 66 , 70], authors adopted the SHO wave function 
to describe a spacial wave function of charmed and charmed- 
strange mesons, where the j3 value in the SHO wave function 
is determined by Eq. (24). This treatment is an approximation 
for simplifying the calculation since the SHO wave function 
'^nLMi(P) is slightly different from the wave function ^(p) in 
Eq. (24). In this work, for the charmed-strange meson A in¬ 
volved in our calculation, we use the numerical wave func¬ 
tion directly obtained as an eigenfunction of the modified GI 
model, which can avoid the uncertainty from the treatment in 
Eq. (24). Eor the mesons B and C in the discussed process, we 
still adopt the SHO wave function, where the corresponding /? 
was calculated in Ref. [36]. 

In the following concrete calculation, the mass is taken 
from PDG [1] if there exists the corresponding experimen¬ 
tal observation of a charmed-strange meson. If the discussed 
charmed-strange meson is still missing, we adopt the theoret¬ 
ical prediction from the modified GI model (the results listed 
in Eig. 2) as the input. 


B. Phenomenological analysis of strong decay 

In this subsection, we caiTy out the phenomenological anal¬ 
ysis of the strong decays of charmed-strange mesons by com¬ 
bining our theoretical results with the corresponding experi¬ 
mental data. 
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TABLE IV: Partial widths of £>*^(2317) —> and Z7ji(2460) —> (in the unit of keV). Here, we list the width of Dji(2460) —> OjTr® 

with a mixing angle dip = -54.7° = - arcsin( yflJS) in the heavy quark limit. 

this work Ref. [90] Ref.[10] Ref. [92] Ref. [93] Ref. [94] Ref. [95] Ref. [44] Ref. [96] Ref. [14] 

r(D*^(23n) ^ D.tt) 11.7 32 21.5 16 ~10 34-44 ^ 6 10-100 155 ±70 3.8 

r(D,,i(2460)^ £>» 11.9 35 21.5 32 ~10 35-51 ^6 - 155±70 3.9 


1. IP slates 

If £)*q( 2317) and Dji(2460) are the IP states in the 
charmed-strange meson family, the OZI-allowed strong de¬ 
cays are forbidden since their masses are below the DKID*K 
thresholds. However, D*g(2317) —» and D5i(2460) —> 

D*7r‘^ can occur via the rj -ti mixing, which are shown in Fig. 
3 [90]. 



FIG. 3: The D*^{2311) 0,71° (left) and £»ji(2460) ^ (right) 

decays through the ij -71° mixing. 

For calculating these decays, we need the Lagrangian in¬ 
cluding the rj - 71^ mixing, which can be expressed as [91] 

= A I- (35) 

4(m„ +md} 

where is a light quark mass matrix, ^ = exp(t7f/f;f) and tt is 
a light meson octet. 

Taking D*q(2317) ^ 0,7^ as an example, we illustrate the 
concrete calculation. The decay amplitude of D*q(2317) ^ 
0,71°^ reads as 

M-'\D*J23\1) ^ D,7i^) = M^Hd*,q( 23\1) ^ D,Ti) 

(36) 

where decay amplitude A1 '^^(Z)*q(2317) ^ 0 , 1 ]) can be cal¬ 
culated by the QPC model. The isospin-violating factor 5,^^! 
is due to the rj - ti mixing, i.e., [97] 

nifi — niu 1 

= - - - - - = -. (37) 

^ m, - (niu + md)l2 43.7 

Adopting the similar treatment, we can derive the formula of 
Dji( 2460) —> Here, Dji(2460) with - 1^ is the mix¬ 

ture between 1 and l^Fi states, which satisfies the relation 


' |D,i(2460)) ' 


/ \ 
cos flip sin flip 


^ |£>.i(2536)) ^ 


-sin flip cos flip ^ 



where the mixing angle dip — -54.7° is obtained in the heavy 
quark limit [36, 98,99]. 


The obtained partial widths of D*q(2317) —> 0,71°^ and 
Dsi(2460) —> 0*71^^ are listed in Table IV. We also compare 
our results with other theoretical predictions from different 
groups [10, 14, 44, 90, 92-96]. Especially, we notice that our 
results are consistent with the corresponding values given in 
Ref [93]. The equality of these two decay widths predicted in 
Ref [14] in the heavy quark limit is well satisfied in our case, 
too. 

Dji( 2536) with - 1^ is the orthogonal partner of 
Dji(2460) as shown in Eq. (38). The D*K channel is its only 
OZI allowed decay mode. In Eig. 4, we present the depen¬ 
dence of the partial decay width of Dji(2536) —» D*K on 
6>ip, which covers the typical value Oip = -54.7° given in the 
heavy quark limit [36, 98, 99]. We notice that our result is 
consistent with the experimental data F = 0.92 ± 0.05 MeV 
[33] and F = 0.75 + 0.23 MeV [34], when the mixing angle 
is around Oip - -52.3° and dip - -53.1°, respectively, which 
are close to dip = -54.7° of the heavy quark limit. In the 
heavy quark limif the decay of Dji(2536) into D*K can occur 
only via a D wave because of the conservation of light degree 
of freedom and hence the decay width is expected to be small 
as shown in Fig. 4. That is, because the experimental mixing 
angle is very close to the one of the heavy quark limit, it is a 
model independent result that Dji(2536) and D*K almost ex¬ 
actly decouples with each other, which has been discussed in 
Refs. [36, 100] 



FIG. 4: The mixing angle dip dependence of a decay width 
Z),i(2536) ^ D*K. 

As a good candidate of the charmed-strange meson with 
1^P2, 0,2(2.513) decays into OK, 0*K, and 0,1] with the par¬ 
tial decay widths 5.42 MeV, 0.57 MeV, and 0.04 MeV, respec¬ 
tively. The total decay width of 0,2(2513) can reach up to 
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TABLE V: Calculated partial decay widths of the OZl-allowed strong decays of 25, ID, 35, and 2P charmed-strange mesons. Units are 
in MeV. Here, forbidden decay channels are marked by - and channels marked by □ are the OZl-allowed modes, which are discussed in 
consideration of the mixing angle dependence. 


Channels £)T2'5o) £)]i(2700) H], 

(2860) 

DAiD(2-)) , 

£)• (2860) Ddi^Po) L»T2P(1 + 

DAW’(2-)) 

)) 0,(2P'(U)) 

OT 2 V 2 ) 

OT3‘5o) 

£> 7335 ,) 

DK - n 

□ 

- 

7.46 

60.86 

- 

- 

2.99 X 10^^ 

- 

23.30 

D*K 76.06 □ 

□ 

□ 

5.98 

- 

62.43 

20.72 

5.30 

38.11 

36.54 

D,r] - □ 

□ 

- 

0.15 

2.07 

- 

- 

0.02 

- 

2.19 

D,t]' 

- 

- 

- 

0.08 

- 

- 

0.001 

- 

0.40 

D*r] - □ 

□ 

□ 

0.06 

- 

1.61 

1.02 

0.18 

2.56 

2.72 


- 

- 

- 

- 

- 

- 

- 

- 

0.07 

DK* 

□ 

□ 

0.37 

- 

39.32 

6.50 

13.46 

3.33 

14.92 

D*K* 

- 

- 

- 

63.66 

48.97 

39.60 

48.32 

18.54 

3.30 X 10-5 

D,(p 

- 

- 

- 

- 

1.45 

8.05 

0.003 

0.03 

0.07 

D*<p 

- 

- 

- 

- 

- 

- 

- 

- 

9.97 X lO-'* 

£)*(2400)7t 

- 

□ 

- 

- 

14.42 

31.13 

- 

38.14 

- 

D*J2-ill)ii - 

- 

□ 

- 

- 

1.06 

2.41 

- 

5.38 

- 

Oi(2430)A: 

- 

- 

- 

3.32 

14.31 

8.11 

14.31 

- 

17.02 

Di{2A2Q)K 

- 

- 

- 

36.16 

19.43 

10.43 

0.75 

- 

5.50 

(2460)7? 

- 

- 

- 

- 

0.25 

0.16 

- 

- 

3.43 

(2536)7? 

- 

- 

- 

- 

- 

- 

- 

- 

0.02 

0^(2460) a: 

- 

- 

- 

- 

82.85 

3.15 

3.90 

5.89 

9.16 

D*J2573)ri - 

- 

- 

- 

- 

- 

- 

- 

- 

0.002 

Total 76.06 

- 

- 

14.02 

166.15 

285.83 

131.28 

86.25 

111.98 

115.35 


6.03 MeV, which is comparable with the experimental value 
16+^ ± 3 [37], 10.4 + 8.3 + 3.0 MeV [41] and 12.1 + 4.5 + 1.6 
MeV [38] given by the ARGUS collaboration and LHCb col¬ 
laboration, respectively. In addition, the branching ratio 


S(D,2(2573) ^ D*°K+) 

—- - - - < 0 33 

S(D,2(2573) ^ D0/:+) 


(39) 


is also measured by the CLEO collaboration [37]. In this 
work, we obtain - 0.106 consistent with 

the present experimental measurement [37]. This value is also 
close to the one, 0.076, predicted in Ref. [14]. 


2. 25 and ID states 

Since there is no candidate for 2'5o charmed-strange me¬ 
son observed in experiment, we adopt the theoretical value 
calculatded by the modified GI model as an input, i.e., 2646 
MeV. As shown in Table V, Dj(2*5o) only decays into D*K, 
where the decay width of Z)j(2'5o) ^ D*K is 76.06 MeV. 

In the following, we shall discuss two observed D*j(2700) 
and D*j(2860) to be an admixture of i) and Ds(l^Di), 

which satisfies the relation 


' |£>:i(2700))' 


cos dsD sin dsD ' 

' 12^51>' 

, |£>:i(2860)) ^ 


^-sin6>5D cos OsD , 



where denotes a mixing angle describing an admixture of 

D,(235i)andD,(l3Di). 

The OsD dependence of the total and partial decay widths of 
D*j(2700) is shown in the Fig. 5. Its main decay modes are 
DK and D*K, both of which were observed in experiments. 
To compare our results with the experimental data, we take 
the BaBar’s measurement in Ref [50] since the width T = 
(149 + together with the ratio listed in Eq. (4) 

were given [50]. As shown in Fig. 5, there exists the common 
OsD range, 6.8° - 11.2°, in which both of the calculated width 
and ratio of D*j(2700) [50] can overlap with the experimental 
data. This small 6s d value is consistent with the estimate of 
OsD in Ref [2]. 

As for the D*j(2860) with 7^ = U recently observed by the 
LHCb collaboration [62, 63] which is considered to be a part¬ 
ner of D*j(2700), the 6s d dependence of the decay behavior 
is depicted in Fig. 6. If taking 6.8° - 11.2° for the range of 
6sD obtained in the study of D*j(2700), we find that the ob¬ 
tained total decay width of D*j(2860) can reach up to ~ 300 
MeV which is comparable with the LHCb data [62, 63] and 
the ratio is S(D;j(2860) ^ D*/:)/S(D;j(2860) ^ DK) = 
0.6 ~ 0.8, which can be tested in future experiment. Our 
study also shows that the main decay channels of D*j(2860) 
are DK{~\AQ MeV), D*K(~95 MeV) and DK*(~50 MeV). 

As an admixture of 1'D2 and 1^D2 states, the states 1D(2") 
and 1D'(2 ) in the charmed-strange meson family satisfy the 













12 





FIG. 7: The 0io dependence of the partial decay widths and total de¬ 
cay width of £)j(l£)(2“)) (left) and £)j(l£)'(2“)) (right). The vertical 
dash lines correspond to the mixing angle 0io = -50.8° 


FIG. 5: The 9sd dependence of the calculated partial and total decay 
widths and ratio T{D*K)IT(DK). The vertical band corresponds to 
the common range of Ogo, where our results can be matched with the 
experimental widths and ratio [50]. 



FIG. 6: The Ogo dependence of the decay widths and ratio 

TiD"K)IT{DK) of D*j(2860). The vertical band corresponds to the 
common range of Qgo, where our results can be matched with the 
experimental width and ratio of D*j(2700) [50]. 


Dj(lD(2 )) and Dj(lD'(2 )) are D*K and DK* . In the heavy 
quark limit, Dj(lD(2 )) couples with D*K via a P wave, while 
it couples with DK* only via a F wave due to the conservation 
of light degree of freedom. As for Ds{lD'{2~)), the situation 
is opposite to D^( 1 D (2 )), Dj(lD'( 2 ')) strongly couples with 
DK* via a P wave. Our present calculations are consistent 
with the conclusion of the heavy quark limit [ 101 , 102 ] and the 
results are model independent in the heavy quark limit. In ad¬ 
dition, The total decay widths of Ds(lD{2~)) and Dj(lD'(2 )) 
can reach up to 240 MeV and 147 MeV, respectively. The 
above study also shows that the D*K and DK* modes are 
the key channels when distinguishing the Dj(lD( 2 ")) and 
Ds{lD'{2^)) states. 

Besides D*^ (2860), the LHCb collaboration also observed 
D*^(2S60). In this work, we test whether 0*^(2860) can be 
a good candidate of Ds(l^Di). The partial and total decay 
widths of D*3(2860) as Dsil^D^) are listed in Table V, in 
which we find that our calculated total width is about one- 
fourth of experimental data [62, 63]. The main decay modes 
of Dsd^D^) are DK and D*K, where the ratio 


SjDAl^D,) ^ D*^K) ^ 

S(D,(dDd D^K) 


(42) 


is obtained. We suggest to do a more precise measurement 
of the resonance parameters and the ratio S{D*'^K)!S{D^K) 
of £>* 3 ( 2860 ), which will finally make a definite conclusion of 
whether D* 3 ( 2860 ) is a Ds(dDd state. 


following relation 


' |l£>( 2 -)> ' 


COS 01 D single ' 


,\ID'{2-)) , 


sin 010 cos 0 io^ 



where 6\d is a mixing angle. In the heavy quark limit, we 
can fix the mixing angle 9\d - -50.8° = - arcsin( V3/5) 
[36, 75, 98]. Adopting the theoretical prediction of the masses 
of Ds{lD{2r)) and Dj(lD'(2 )) as input, we list their allowed 
decay channels in Table V. The 9\d dependence of the par¬ 
tial and total decay widths of Dj(lD(2“)) and D^(1D'(2 )) is 
given in Fig. 7. When taking the limit value 9xo - -50.8° 
[36, 75, 98], we conclude that the main decay modes of 


3. 2P stales 

There are four 2P states. Among these, Ds{2?Pq) and 
Ds{ 2 ^P 2 ) are still missing in experiment, while there exist pos¬ 
sible candidates for Ds(2P(d)) and Ds(2P'(d)). 

The partial and total decay widths of Dsd^Po) are shown 
in Table V, where the mass of Ds{2^Po) is fixed as 

2960 MeV (see Table. II). The total decay width 
ofDs(2^Po) is 166.15 MeV and there exist three main decay 
modes, DK, D*K*, and D\{2A2Q)K (see Table. V). 

As for the Dsd^Pi) state with the predicted mass 3004 
MeV, the main decay modes are D*K*, D*K, and Di(2430)£r 
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with the total decay width 86.25 MeV. The ratio is predicted 
as S(DA 2 ^P 2 ) ^ D*K)IS(DA 2 ^P 2 ) ^ DK*) = 0.39. 

We notice that there is an evidence of the structure around 
2960 MeV in the invariant mass spectrum given 

by LHCb [62, 63] except for the observed D*j(2860) and 
£>* 3 ( 2860 ). If this evidence will be conhrmed by future ex¬ 
periment, this structure around 2960 MeV must be either 
D,( 2 ^Po) or D,( 2 ^P 2 ) since there are no other candidates 
for natural states from 2800MeV to 3100 MeV as shown 
in Fig. 2. According to the calculated decay behaviors of 
D,( 2 ^Po) and D,{ 2 ^P 2 ), we can exclude the possibility of the 
Ds{ 2 ^P 2 ) assignment to this structure since the decay width of 
Ds( 2 ^P 2 ) DK is quite small. 

In the following, we discuss the possibilities of the observed 
£>*^(3040) as D,(2£’(l+)) or D,(2£’'(l+)), where D,(2P(1+)) 
or Dj( 2 £’'(l+)) satisfy 


' |2F(l+)> ' 


cos 02 P sin 02 P 

'|2'Pi)' 

, |2F'(1^)> , 


^ - sin 02 P cos 02 P ^ 



with a mixing angle 02 p. If taking 62 P - 0ip - -54.7° [36, 98] 
to estimate the decay behaviors of Dj( 2 £’(l^)) or Dj( 2 £’'(l^)), 
we obtain the numerical results listed in Table V. 

If assigning D*y(3040) to Dj(2P(l^)), the calculated total 
decay width is 285.83 MeV, which is consistent with the ex¬ 
perimental width [50]. The main decay modes are D*K, DK*, 
D*K* and D*,^{2A6Q)K, which can explain why D*y(3040) was 
first observed by experiment in the D*K channel [50]. Addi¬ 
tionally, we predict the ratio 


S(D,(2£’(1+)) ^ D*K) 
S(D,(2£’(1+)) ^ DK*) 


(44) 


As for D 5 ( 3 '5 0 ), the total decay width is about 111.98 MeV 
and the main decay channels are D*K, D*K* and D*^^{2AQQ)K. 
As for Ds(3^S 1 ), the total decay width is 115.35 MeV and 
the main decay channels are DK, D*K, DK* and Di(2430)£r. 
These informations are valuable for experimental search for 
these two missing Dj(3'5o) and DsO^S 1 ) mesons. 

As the unnatural state, the predicted mass (3092 MeV) of 
Dj(3'5o) is very close to the mass of D*^(3040) reported by 
the Babar Collaboration [50]. Thus, we again check whether 
D*^(3040) is explained as Ds{3^So). Since D*^(3040) was ob¬ 
served in the D*K channel which is one of the main decay 
channels of 0^(3*5 0 ), this fact cannot exclude this possibility. 
However, the obtained total decay width of Dj(3*5o) is a lit¬ 
tle bit smaller than the lower limit of the experimental width 
of D*^(3040). Although the predicted Ds(3^So) can not fit 
all the experimental feature of D*^(3040), there exists small 
possibility of D*y(3040) as Dj(3'5o) if considering large ex¬ 
perimental error of D*j(3Q4-Q). Thus, a crucial task in future 
experiment is the precise measurement of the resonance pa¬ 
rameters of D*^(3040), which can help us to have a definite 
conclusion of the properties of D*^(3040). 


5. 2D states 

There are four 2D states all of which are absent in exper¬ 
iment. In Table VI, we predict their decay properties. As 
for the Ds{2^D\) state, the predicted mass by the modihed GI 
model is 3244 MeV and the total width obtained in the QPC 
model is 165.20 MeV. Its main decay channels are DK, D*K 
and D\(2A20)K and the typical decay ratio is given by 


which can be tested by future experiment. The above study 
shows that D*^(3040) is a good candidate of Ds(2P{V)), 
which is the hrst radial excitation of Dji(2460). 

We need to discuss another possible assignment of 
£)*y(3040) to D(2£’'(l^)) if considering study of only the mass 
spectrum shown in Fig. 2 because the mass of D*^(3040) is 
very close to the predicted mass of Ds{2P'{\^)). Moreover, 
the total decay width of D^(2£’'(l+)) is 131.28 MeV, which 
is comparable to the lower limit of experimental width of 
D*y(3040) [50]. The predicted ratio is, however, 

nDA2P'(l-))^D*K) 

B{Ds(2P'(\+)) ^ DK*) 

which is quite different from that given in Eq. (44). 
Thus, we suggest to carry out the measurement of the ratio 
S(D;^(3040) ^ D*K)IB{D*J3QAQ) DK*), which will be 
helpful to hnally identify the inner stmcture of D*^(3040). 


4. 35 states 

In this subsection, we predict the decay behaviors of two 
35 states in the charmed-strange meson family, which are still 
missing in experiment. As shown in Table V, the partial and 
total decay widths of Dj(3'5o) and Ds{3^S 1 ) are calculated. 


S{Ds{2^Di)) ^ D*K) 
S(D,(23Di)) ^ DK) 


(46) 


As an admixture of 2 'D 2 and 2^D2 states, Ds(2D(2 )) and 
Ds{2D'{2~)) satisfy the following relation 


' |2£»(2-)) ' 


cos 02D sin 02D 

' |2'£»2> ' 

, |2£»'(2-)> ^ 


^ - sin 02D cos 02D , 

. |2'£)2> , 


(47) 


where the mixing angle can be fixed as 02 d - -50.8° in the 
heavy quark limit [36, 75, 98] when discussing their decay 
behaviors. 

The Ds(2D(2^)) with the predicted mass 3238 MeV has 
the total width 141.49 MeV and its main decay channels are 
D*K and D*^,^{2A6Q)K that contribute to almost 80% of the to¬ 
tal decay width. As the partner of Ds{2D{2^)), Ds(2D'(2^)) 
has the predicted mass 3260 MeV, is a narrow state com¬ 
pared with Ds{2D(2~)) because the total decay widths is 52.4 
MeV, and its main decay mode is DK* with the partial decay 
width 30.67 MeV. The difference of the total decay widths 
of Ds(2D(2~)) and Ds(2D'(2^)) can be understood since the 
decays of Dj(2D(2')) and Ds{2D'{2~)) into D*K occur via P- 
wave and F-wave in the heavy quark limit [75], respectively. 

As for Dsii^Di,) with mass 3251 MeV, the total decay width 
is 44.85 MeV and the main decay channels are D*K* and 

















14 


D\(2A30)K. Here, the ratio 


The D,(1T’(3+)) and D,(1T'(3+)) satisfy 


S{Ds{2^Di)) ^ D*K) 
S(D,(23D3)) ^ DK) 


(48) 


' |l^(3+)> ' 

/ 

, |1F'(3+))^ 

\ 


cos 0\F 
- sin OiF 


sin 6 \f 
cos OiF 




(50) 


is predicted via the QPC model. 


TABLE VI: Decay behaviors of four 2D charmed-strange mesons (in 
the unit of MeV). 


Channels 

D,{2^Di) 

DA2D(2-)) D,(2D'(2-)) 

DA2?Di) 

DK 

65.02 

- 

- 

5.25 

Dst] 

4.13 

- 

- 

0.23 

D.ri' 

0.87 

- 

- 

0.02 

D'K 

20.32 

76.78 

0.03 

0.96 

Dlrj 

0.93 

4.12 

0.006 

0.04 

D*srj' 

0.04 

0.29 

0.002 

3.0 X 10-’ 

DK* 

5.03 

4.73 

30.67 

0.69 

D,(p 

0.06 

0.12 

0.04 

0.04 

D*K* 

3.2 X 10-* 

1.31 

3.11 

19.64 

d:<p 

0.37 

0.42 

0.47 

0.36 

D*{2400)K 

- 

6.82 

8.00 

- 

D:j23n)ri 

- 

0.36 

0.37 

- 

D*g(2400)K* 

0.54 

0.16 

0.53 

0.63 

Di{2430)K 

14.25 

8.14 

5.64 

11.50 

Dii2420)K 

39.90 

2.80 

1.71 

1.26 

Dsi(2460)rj 

0.59 

0.32 

0.24 

0.43 

D,i{2536)ti 

0.19 

0.07 

0.06 

0.007 

D*(2460)/f 

12.77 

34.22 

1.46 

3.74 

D*;j2513)ri 

0.19 

0.83 

0.06 

0.05 

Total 

165.20 

141.49 

52.40 

44.85 


6. \F states 


Similar to the experimental situation of the 2D states in 
the charmed-strange meson family, four IF charmed-strange 
mesons are still missing in experiment. Thus, in the following 
we predict their decay behaviors. 

The mass of Dj(l^F 2 ) is 3159 MeV predicted by the modi¬ 
fied GI model. The OZI-allowed two-body strong decay chan¬ 
nels are listed in Table VII. Our calculation shows that its 
main decay modes are DK, D*K, DK*, and Di(2420)F. The 
DsO^F 2 ) is a broad state with the total decay width 415.97 
MeV, where the predicted ratio is 


S{DA2^F2)) ^ D*K) 
S{DA2^F2)) ^ DK) 


(49) 


where 6 if is the mixing angle, which can be determined as 
OiF - -49.1° = - arcsin(2/V7) in the heavy quark limit 
[36, 98]. In Table VII, we collect the calculated decay widths 
of Z)5(1F(3^)) and Dj(lF'(3^)), where we take 3139 MeV 
and 3169 MeV as their mass input, respectively. The main de¬ 
cay channels of D,(1F(3+)) are D*K, D*K* and, D*(2460)F 
and the total decay width can reach up to 372.48 MeV. Thus, 
Dj(lF(3^)) has a very broad width. As for Di(lF'(3^)), DK* 
and D*K* are its main decay modes and the total width of 
Dj(lF'(3+)) is 193.34 MeV. 

The Dsil^F 4 ) with the mass 3143 MeV dominantly decays 
into D*K* and the total decay width of Dj(PF 4 ) is 150.79 
MeV. In addition, we predict the following ratios 


S{DA2^F4)) ^ DK) _ 
B{D,0F4)) ^ D*K) 

and 

S{Dsi2^F4)) ^ DK*) ^ 
S(D,(23F4)) ^ DK) 
which can be tested in future experiment. 


(51) 


(52) 


V. SUMMARY 

In the paste decade, more and more charmed-strange 
states have been reported by different experiments, which 
are collected in Table I. The present status of these ob¬ 
served charmed-strange states has stimulated us to revisit the 
charmed-strange meson family and to systemically carry out 
the study of their mass spectra and two-body OZI-allowed 
strong decays. 

In this work, we have adopted the modified GI model to get 
the mass spectra of charmed-strange meson family, where the 
screening effect partly reflecting the unquenched effect is con¬ 
sidered in our calculation. Comparing the theoretical results 
with the experimental data, we can roughly obtain the prop¬ 
erties of the observed charmed-strange mesons. In addition, 
we have also predicted the masses of the higher radial and or¬ 
bital excitations in the charmed-strange meson family. This 
information is important in searching for these missing higher 
charmed-strange mesons in future experiment. Besides the 
mass information, in our calculation we have obtained the nu¬ 
merical results of the spatial wave functions of the discussed 
charmed-strange mesons, which are applied to calculate their 
two-body OZI-allowed strong decays. 

To obtain the decay behaviors of the discussed charmed- 
strange mesons, we have adopted the QPC model in this work, 
where the IF, 2P, ID, 25, 2D, 35, and IF states in the 
charmed-strange meson family are involved. The study of 
their strong decay behaviors further test the possible assign¬ 
ments to the observed states, where a phenomenological anal¬ 
ysis is performed. As for the higher charmed-strange mesons 
















15 


TABLE VII: Decay behaviors of four IF charmed-strange mesons ( 
in the unit of MeV). 


Channels 


DT1C(3+)) 

DT1C'(3G) 


DK 

57.82 

- 

- 

5.30 

Dst] 

3.16 

- 

- 

0.11 

Dstl' 

0.41 

- 

- 

0.002 

D*K 

43.12 

98.10 

14.86 

6.88 

D'sV 

2.19 

4.98 

0.28 

0.11 

DW 

0.06 

0.08 

4.3 X lO-'* 

5.5 X 10-^ 

DK* 

41.12 

9.34 

102.32 

4.35 


0.85 

0.04 

2.36 

0.08 

D*K* 

16.68 

51.12 

66.16 

130.87 

D'A 

0.01 

0.002 

0.09 

0.01 

D5(2400)K 

- 

0.50 

3.32 

- 

£>:o(2317)7, 

- 

0.02 

0.16 

- 

Di(2430)K 

0.49 

0.22 

0.68 

1.24 

Di(2420)K 

220.10 

1.40 

2.14 

0.42 

D^i (2460)7/ 

0.007 

0.003 

0.02 

0.02 

D,i(2536)ri 

5.35 

6.8 X 10-“* 

0.004 

3.0 X 10-^ 

D*(2460)K 

24.34 

205.68 

0.95 

1.40 

0^2(2573)7/ 

0.26 

0.99 

7.9 X lO-'* 

6.5 X 10-5 

Total 

415.97 

372.48 

193.34 

150.79 


absent in experiment, we have predicted their partial and total 
decay widths, and some typical decay ratios, which is valuable 
for experimental study of these states. 

Since 2003, the BaBar, Belle, CLEOc, and LHCb experi¬ 
ments have made a big progress on the search for charmed- 
strange mesons. In the next ten years, we believe that more 
candidates of the charmed-strange mesons will be reported 
with running of the LHC experiment at 14 TeV collision en¬ 
ergy, forthcoming Bellell and PANDA. The study presented 
in this work is helpful to identify these observed charmed- 
strange states and carry out a search for higher radial and or¬ 
bital excitations in the charmed-strange meson family. 
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